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Abstract 

We explore net charge probability distributions in heavy ion collisions within the 
hadron resonance gas model. The distributions for strangeness, electric charge and 
baryon number are derived. We show that, within this model, net charge probability 
distributions and the resulting fluctuations can be computed directly from the mea- 
sured yields of charged and multi-charged hadrons. The influence of multi-charged 
particles and quantum statistics on the shape of the distribution is examined. We 
discuss the properties of the net proton distribution along the chemical freeze-out 
line. The model results presented here can be compared with data at RHIC energies 
and at the LHC to possibly search for the relation between chemical freeze-out and 
QCD cross-over lines in heavy ion collisions. 
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1 Introduction 



Detailed studies of hadron production in nucleus-nucleus collisions from SIS 
up to LHC energies provide the very interesting result that particle yields 
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exhibit thermal characteristics and are well described by the statistical oper- 
ator of a hadron resonance gas (HRG) [Tp|3] . Furthermore, there is a unique 
relation between the collision energy and the corresponding thermal param- 
eters, which defines the so called chemical freeze-out line in the temperature 
and chemical potential (T, /i)-plane [T|2|H[5] . The phenomenological freeze-out 
curve, along which particle yields attain their measured values, is considered 
to be indicative for the QCD phase boundary. 

Moreover, results of lattice QCD studies on the critical temperature at 
vanishing chemical potential piT] and on the slope parameter characterizing 
the dependence of Tc on /i pi9] show that, at high energies, the chemical 
freeze-out curve is very close to the QCD (cross-over) phase transition. Al- 
though the origin of this fact is not fully understood, the observation is sup- 
ported by dynamical arguments [TU] linking the phase transition to chemical 
freeze-out. The HRG partition function also describes the thermodynamics of 
a strongly interacting medium in the hadronic phase as obtained in lattice 
QCD [TTllTMMllTSllTGlITT] . 




The proximity of the freeze-out curve to the phase boundary suggests that 
the QCD phase transition and its related critical properties should be ob- 
servable in heavy ion collisions. However, an experimental verification of a 
phase change in a medium created in such collisions requires sensitive probes. 
In this context, a particular role is attributed to fiuctuations of conserved 

Recently, it was argued that, at high energies, the history of the system, in 
particular the path through the QCD cross-over transition from the deconfined 
and chirally symmetric phase to the hadronic one, may be refiected in fiuctua- 
tions of conserved charges, specifically in their higher cumulants [2^24|25|I2B] . 
The characteristic signature of such transition may thus be manifested in de- 
viations of higher cumulants of the charge distributions from the HRG results, 
if the freeze-out happens near or at the QCD phase boundary. At vanishing 
chemical potential, the sixth and higher order cumulants can be negative, also 
in the hadronic phase, in contrast to predictions from HRG model, where only 
positive yields are obtained [23]. It is therefore useful to consider the HRG 
model results on moments of charge fiuctuations as a theoretical baseline; any 
deviation from this could be an indication for critical phenomena at the time 
of hadronization [TT1I231I21] . 

First data on cumulants of the net proton multiplicity were recently obtained 
by the STAR Collaboration in Au-Au collisions at several energies and central- 
ities [27p8] . The basic properties of the measured fiuctuations, in particular 
of their ratios, are consistent with HRG model expectations [23|28] . 




In practice, the cumulants of the net proton fiuctuations are obtained from the 
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corresponding probability distribution [3T|32f33|IM] . As recently noted [35], a 
direct comparison of the net proton multiplicity distribution, measured in 
heavy ion collisions, with that of the hadron resonance gas may provide a 
deeper insight into the relation between chemical freeze-out and the QCD 
cross-over transition. The comparison is facilitated by the direct link between 
the net proton probability distribution and the mean number of protons and 
antiprotons in the HRG model |35]. This eliminates uncertainties connected 
with the extraction of freeze-out conditions and thus provides a less ambiguous 
method for confronting the model results with data. 



A freeze-out close to the phase boundary should be reflected in a modification 
of the width and shape of the net proton probability distribution [33]. In the 
crossover regime from hadronic matter to the quark gluon plasma the width 
of the distribution is expected to be reduced, compared to the HRG model 
predictions. If, on the other hand, freeze-out occurs near the chiral critical end 
point, the width could be enhanced relative to that obtained in the hadron res- 
onance gas. The first STAR results for the net proton probability distribution 
from Au-Au collision at i/ijviv ~ 200 GeV suggest that the width of the ex- 
perimental distribution may indeed be reduced [35] . In order to scrutinize this 
effect, a detailed analysis of data on the net-proton probability distributions 
at different beam energies is necessary. 



In this paper we extend our previous results [35] and examine further the 
significance of the net charge probability distributions as a tool for exploring 
critical phenomena in heavy ion experiments. In particular, we consider the 
probability distributions related to the fluctuations of strangeness and electric 
charge. We show that, in the hadron resonance gas, these distributions, as 
well as all corresponding cumulants, can be directly linked to mean particle 
multiplicities. The influence of multi-charged hadrons and quantum statistics 
on the shape of net charge distributions is discussed. We compute the net 
proton distribution along the freeze-out line and study its systematics. 



In the next section we derive probability distributions of conserved charges 
in a thermodynamic system and apply the results to the hadron resonance 
gas. In Section III we discuss the properties of these distributions in the HRG 
model and make predictions for the net proton probability distribution and 
the corresponding fluctuations in heavy ion collisions at RHIC energies and 
at the LHC. Finally, in Section IV we give our conclusions. 
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2 Probability distribution of conserved charges 



2.1 General case 



In statistical physics the probability distribution P{N) can be obtained di- 
rectly from the relation between the canonical Z{N, T, V) and grand canonical 
E{lJ,q,T,V) partition functions. The thermodynamics of a system with fixed 
net charge N is described by Z{N,T, V), which is obtained from the density 
operator 

Z{N,T,V)=TTMe-^". (1) 
Here the trace is constrained to configurations with a given net charge A^. The 
grand canonical partition function 

Z(/i„r,V) = Tre-^^+'^^, (2) 

where fi = f^q/T, is related to Z{N,T, V) by the cluster decomposition 

Z{fiq,T,V)=J2Z{N,T,V)ef'''. (3) 

N 

The physical meaning of each term in the above sum is the probability of 
finding a configuration with net charge N. With the proper normalization, 
the probability is [T|31f32j 

P{N) = \,,, ZiN, T, V)e^''. (4) 

By eliminating Z{fiq,T,V) in favor of the thermodynamic pressure InZ = 
VT^p{T, fi), one arrives at the probability distribution 

P{N) = Z{N, T, y)eA^-^^'p(T,A) 

for finding a state with net charge in a system of volume V at the tempe- 
rature T. 

The canonical partition function can be directly obtained from the grand 
canonical one computed at imaginary chemical potential [1], 

Z(N, T,V) = — r dde^'^^ZiiTd, T, V), (6) 
27r JO 

where the chemical potential was Wick rotated by the substitution fi — > iO. 
Equations ^ and ([6]) define the probability distribution of net charge in a 
thermal system of volume V . 

The n-th order cumulant Xn of the net charge is obtained by computing 
the corresponding moments of P{N). Alternatively, the cumulants are given 
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Fig. 1. Left-hand and middle figures: The pf-distributions of antiprotons and pro- 
tons in Au-Au cohisions at -v/i^rjv =130 GeV for different centrahties [36] ■ The hnes 
are blast-wave fits [37] to the data. Right-hand figure: The net-proton distribution 
in the hadron resonance gas computed from Eq. (jlSp for different centralities. The 
full lines are obtained using experimental inputs for the proton Np and antiproton 
Np yields, while the broken-lines are obtained with Np and Np computed in the 
thermal model with parameters describing the chemical freeze-out |36j . 

by the generalized susceptibilities 

Thus, we have for the susceptibihty or the second order cumulant, 

where the moments of P{N) are defined by 

{N'') = Y,N''P{N). (9) 

N 



2.2 The net-charge distribution in a hadron resonance gas 



In a strongly interacting medium there are conservation laws related to the 
global symmetries of the QCD Lagrangian. In the following, we explore prob- 
ability distributions of such conserved charges and the corresponding fluctua- 
tions. We focus mainly on fluctuations of the net baryon number in hot and 
dense hadronic matter, but give results also for fluctuations of strangeness and 
electric charge. 

We use the hadron resonance gas partition function to model the thermody- 
namics of strongly interacting matter. This model yields a good approximation 
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to the lattice QCD equation of state below the QCD phase transition. In the 
HRG model the thermodynamic pressure is a sum of meson and baryon contri- 
butions. The conservation of baryon number, strangeness and electric charge 
is accounted for by introducing the corresponding chemical potentials. 

In Boltzmann approximation, we can write the thermodjTiamic pressure in 
the following form: 

kl 

PVp{T,fi) = YX^ne""^ + ^ne-"'^]. (10) 

n=0 

Here one of the conservation laws is explicitly exhibited. We denote the corre- 
sponding chemical potential generically by fx. Depending on which conserva- 
tion law is picked, baryon number, electric charge and strangeness, the maxi- 
mum charge q take the values |g| = 1, 2 or 3, respectively. The remaining con- 
servation laws are implicitly included in the parameters Zn and which en- 
code the thermal phase-space of all particles and antiparticles carrying charge 
n and — n, respectively. Within the HRG model, we find 

VT 

zn= E 4> 4 = ^^?.^'^2K/T)e^'-'^^ (11) 

iShadrons 

where the sum is taken over all stable hadrons as well as all hadron resonances. 
Moreover, the vector /j represents those charges carried by particle i that are 
not explicitly shown in Eq. ( ITOll . while /Tg subsumes the corresponding chemical 
potentials. The Zn is obtained from Eq. (fTTl) by the replacement /2 — )■ — /i. 

The phase-space parameters Zn and related to the mean-number of 

particles Nn and antiparticles iV_„ carrying charge n and — n, respectively, 

iV„ = 2;„exp(n/i), iV_n = ^„exp(-n/i). (12) 

The partition function, where one conservation law is treated exactly, is ob- 
tained using Eqs. (jH]) and ffTU]) . We denote the resulting partition function 
by where |g| is the relevant maximal charge. We start with the most 

complex case, i.e. with |g| = 3. Within the HRG model the corresponding 
canonical partition function is given by [1] 

oo oo 

Z^^\N,T,V)= J2 E ar'^"'''"'^4«3^^-2.-3fc(xi)/.(x2)/,(x3), (13) 

i=—oo k=—oo 

where In{x) is the modified Bessel function. In Eq. f|T3l) . we have also intro- 
duced the notation 
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On = J — , and Xn = 2v^i^ = 2 y'AyvI^ (14) 
y 

with n G (1, 2, 3). Within the hadron resonance gas model the canonical par- 
tition function Z^^^ is appropriate for describing a system with exact conser- 
vation of strangeness. 

For the electric charge, the appropriate canonical partition function is Z^'^\ 
since there are single- and double-charged hadrons. We stress, however, that 
this is correct only within the Boltzmann approximation we consider here. In 
the electric charge sector this approximation is not suitable, as we will discuss 
later. The generic form of Z^'^^ is obtained from Eq. f lT3|l by taking the limit, 
X3 — 7- and 03—7-1. This yields a non-zero contribution only for k = 0, with 
lima.3_>.o /A:=o(a;3) = 1, leaving the summation over i. The canonical partition 
function for baryon number conservation, Z^^-*, is then obtained from Z^^-* by 
taking the corresponding limit, X2 ^ and 02 — >■ 1. Again, only the i = 
contribution survives, leaving an expression in closed form [31] 

\N/2 

zW(iV,T,y) = (^) lN{2^n). (15) 



In the HRG model, the probability P(N) to find a state with net-charge N is 
obtained from Eqs.(jl]), ffTOj) and f|T3|l . We compute explicitly the probability 
distributions relevant for the phenomenology of heavy-ion collisions, i.e. for 
strangeness (5"), electric charge {Q) and baryon number (5). 

Since there are hadrons with strangeness one, two and three, the relevant 
probability distribution is that obtained from Z^'^^: 



P{S) 



_ , s 



exp 



.n=l 



i/2 



'A' 



/. 2JS0S 



00 00 

E E 

-i-3fc/2 



■^3 



fc/2 



h 2x S^S. 



'3^-3 



W-5 2JSiS^i , (16) 



where Sg and S-s denote the mean-number of particles and antiparticles car- 
rying strangeness s, where s G (1,2,3) (see Eq. (|T2l) ). 

For the electric charge, treated in Boltzmann approximation, the correspond- 
ing probability is given by 
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PiQ) 



exp 



.n=l 



E 



'Q2 



i/2 



1 ' 



(17) 



where Qg and Q^g denote the mean number of particles and antiparticles with 
electric charge q and —q respectively, with q G (1,2). 

Finally, the probability P{B) of net-baryon number B can be expressed in 
terms of the mean-number of baryons Bi and antibaryons [35] by the 
Skellam distribution. 



P{B) 



'A' 

B-i 



B_ 

2 



Ib 2JBiB^i exp Bi + B_i 



(18) 



Equations (fT6!) - (fT8ll determine the probability distributions for the fluctu- 
ations of net strangeness, electric charge and baryon number in a hadron 
resonance gas in equilibrium and under Boltzmann statistics. We note that 
the probability distributions of the HRG depend only on the mean-number of 
charged particles and antiparticles. 



3 Net-proton distribution in heavy ion collisions 



3.1 Probability distributions 



The mean number of particles Nn and antiparticles N^n obtained in experi- 
ment can, when available, be used directly in Eqs. ( IT614T81) . where the volume 
and all thermal parameters have been eliminated. Thus, predictions of the 
HRG model can be compared with data in an unambiguous way, avoiding fur- 
ther model assumptions. In particular, using measured multiplicities Nn and 
N-n has the advantage that experimental cuts are included consistently for all 
quantities. We note that for a comparison of data to theory the experimental 
results have to be corrected for effects due to overall baryon, strangeness, and 
charge conservation. 

Alternatively, in the HRG model the mean number of particles and antipar- 
ticles N_n entering P{N) in Eqs. (fTGHTSl) can, given the thermal and volume 
parameters, be computed using Eq. (|T2i) . In applications of the HRG model to 
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heavy-ion phenomenology the thermal parameters are determined along the 
chemical freeze-out curve. The volume can then be fixed to reproduce the net 
charge number or particular particle yields. 

At present, the first approach can be applied only for the net baryon multiplic- 
ity distribution. The electric charge and strangeness probability distributions 
require experimental input on yields of all charge states in the acceptance win- 
dow where the probability distributions are obtained. At present, such data 
are not available for strangeness and electric charge. On the other hand, first 
data on charge fiuctuations and higher order cumulants of net proton multi- 
plicities were recently obtained by the STAR Collaboration in Au-Au collisions 
at several collision energies [27p8] . The data were taken at mid-rapidity in a 
restricted range of transverse momenta, 0.4 GeV < < 0.8 GeV. 

Consequently, in the following, we focus on the net-proton probability distri- 
butions in the hadron resonance gas model. 

The probability distribution, Eq. f|T8|) . is readily generalized to protons in a 
momentum window. Due to the factorization of the grand canonical single- 
particle partition function in momentum space, Eq. yields a normalized 
probability distribution, when both partition functions are restricted to the 
same momentum window. Similarly, the limitation to protons follows from the 
factorization of the contribution of any particle species to Z. Consequently, 
in the HRG model the probability distribution of net proton number in a 
momentum window, 6Np = Np — Np, is obtained from Eq. fll8l) by simply 
replacing the mean number of baryons Bi and anti-baryons by that of 
protons Np and antiprotons Np in the same kinematic window. 

Thus, we can compute the net-proton distribution using Eq. f fTSl) . provided we 
have access to the mean values Np and Np measured in the same kinematic 
window. Recently, such studies were performed in Ref. [35], using the STAR 
data [271136] on prdistributions of antiprotons, the mean number of net protons 
as well as P{N) for Au-Au collisions at ^/sNN = 200 GeV and at several 
centralities. To assess the dependence of P{N) in A-A collisions on beam 
energy, we extend these studies to different energies along the chemical freeze- 
out curve. 

Unambiguous predictions of the model for the net proton probability distri- 
bution can be obtained in Au-Au collisions at y/sj^j^ = 130 GeV as well as 
in Pb-Pb collisions at the LHC, since at these energies the pj-distributions of 
protons and antiprotons are available [5B|55] . 

In Fig. [H we show the prspectra of protons and antiprotons at y/sj^j^ = 130 
GeV obtained by the STAR Collaboration [36] • By fitting these spectra with 
the blast- wave model [37] and integrating the data in the range 0.4 GeV < 
Pt < 0.8 GeV, we find Np and Np for the different centralities. 
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Fig. 2. Left-hand figure: The pj-distribution of antiprotons obtained by the STAR 
Collaboration in Au-Au collisions at \/s7V7V ~ ^2.4 GeV for different centralities 
|36| . The lines are blast- wave fits [S^ to the data. Right-hand figure: The net-pro- 
ton distribution calculated in the hadron resonance gas from Eq. (jl8|) for Au-Au 
collisions at \fsj\fj^ = 62.4 GeV at different centralities. The full-lines are obtained 
with experimental inputs for the proton Np yields and the mean net baryon number 
M = Np — Np, while the broken-lines are obtained with Np and Np calculated in 
the thermal model with parameters describing the chemical freeze-out |36j . 

In the right panel of Fig. [1] we show the corresponding HRG result for the 
net proton multiplicity distribution in Au+Au collisions at ^/SNN =130. Also 
shown in this figure are the results for P{N) obtained within the HRG model 
with Np and Np computed at the chemical freeze-out using the thermal pa- 
rameters of Ref. [36] • The volume was fixed by requiring that the measured 
mean-number of protons be reproduced. The overall agreement of the two 
methods is satisfactory. Nevertheless, a direct calculation of P{N) from the 
measured yields of protons and antiprotons is preferable, since it is unaffected 
by systematic errors in the thermal parameters. 

In order to find P{N) in Au-Au collisions at a/s^^ = 62.4 GeV, we used the 
measured prspectra of antiprotons [36] to determine the antiproton multiplic- 
ity. The proton yield was then fixed using the measured mean net proton num- 
ber [27|28] . We stress that the data on net-proton number are not corrected 
for efficiency [27f29] . This implies a systematic uncertainty in our results. The 
resulting distributions for different centralities are shown in Fig. |2J Again, we 
also show results of the thermal model, where Np and Np are computed at 
the chemical freeze-out using the thermal parameters of Ref. [36] • The two 
methods give similar results on P{N) at all centralities. 

Recently, the STAR Collaboration has also presented preliminary data on mo- 
ments of the net-proton distribution in central Au-Au collisions at y/sj^j^ =7.7, 
11.5 and 39 GeV. However, at these energies the corresponding p^-spectra of 
protons and antiprotons are not available. Thus, for these energies we com- 
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Fig. 3. Predictions of the hadron resonance gas model for the net proton probabihty 
distribution in central nucleus-nucleus collisions at different energies. The results 
for Au-Au collisions at -v/s^vAf ~ 4 and y/sj\ij^ = 130 GeV are obtained using 
Eq. (jlSp with experimental inputs for the proton Np and antiproton Np yields 
extracted from the pt spectra of the STAR Collaboration [HU]. The probability 
distribution for Pb-Pb collisions at -s/sjviv ~ ^•''''^ '^^^ obtained using Eq. (fT8|) 
with experimental inputs for Np and Np extracted from the preliminary p^-spectra 
of the ALICE Collaboration [38]. At lower energies, P{N) is computed with Np 
and Np determined at the chemical freeze-out with thermal parameters from Ref. 
[2]. The required experimental inputs on mean net-proton numbers are from STAR 
Collaboration |27l28j . 

puted P{N) in the hadron resonance gas using the proton and antiproton 
yields obtained in the thermal model at chemical freeze-out with the thermal 
parameters of Ref. [2]. The volume was fixed by fitting the mean number of 
net protons to the STAR data [30]. The resulting probability distributions 
P{N) for these three energies are shown in Fig. [31 

In order to assess the possible role of dynamical effects on the probability dis- 
tribution, we compare in Fig. IHthe STAR data on P{N) in Au-Au collisions 
at v^TVAf =200 GeV and the HRG results of Ref. [35] with the net proton prob- 
ability distribution obtained with the UrQMD event generator. Deviations of 
the thermal model from the data and their possible origin were discussed in 
Ref. [32] • We note that the distribution obtained in the UrQMD model is 
much narrower than the data. The reason for this remains to be clarified. 

In Fig. [3] we summarize the HRG model results for the net proton distributions 
in central nucleus-nucleus at energies spanning from y/sj^j^ =7.7 GeV up to 
y/sj^j^ =2.74 TeV. The probability distribution at the LHC energy is computed 
with proton and antiproton yields obtained by integrating the preliminary data 
on pf-spectra of the ALICE Collaboration [38] in the pt range (0.4 GeV < 
Pt < 0.8 GeV). The P{N) distributions shown in Figs ([T]l3]) can be compared 
directly to properly normalized data from RHIC and LHC. 
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3.2 Beam energy systematics of P{N) and its cumulants in heavy ion colli- 
sions 



As seen in Fig. |3l there are systematic changes in the net-proton probabihty 
distributions in central heavy ion colhsions with increasing beam energy. There 
is a clear asymmetry in P{N) at lower energies, which is largely eliminated 
at the LHC. With increasing y/sj^j^, there is a corresponding shift of the 
maximum of P{N) from a nonzero value of towards = 0. These properties 
of P{N) are direct consequences of the decreasing net baryon density with 
increasing beam energy. 

A rather striking feature, seen in Fig. |31 is that in the entire energy range 
from the lowest RHIC energies up to LHC energies the maximum of P{N), 
is almost independent of energy. This is unexpected, as it is well known that 
in heavy ion collisions the total density of baryons plus antibaryons in the 
whole pt-range at mid-rapidity is approximately constant along the chemical 
freeze-out curve [32] • The yield of both protons and anti-protons thus changes 
considerably with energy. On the other hand, as we show in the following, the 

maximum of P{N) is proportional to Np + Np 

The maximum value of P{N) can be obtained analytically from Eq. (fT8|) in 
two limiting cases: for very high and very low energies. The probability P{N) 
has its maximum at some value of A^ ~ [(A^)], where [A^] denotes an integer 
part of A^. 

At high energy, the mean net-proton number is small, (A^) = 0, while Np is 
large. Therefore, the maximum of P{N) appears at A^ ~ [(A^)] = and 



-{Np+Np) 



271 {Np + Np 



-1/2 



(19) 



where we used the asymptotic expansion of the Bessel function 



IJx u + 1) 



V2 



Tlx 



(20) 



as well as that Np ^ Np and 2 



Np ■ Np 



1/2 



{Np + Np 



At low energies, on the other hand, the number of antiprotons is very small 
Np <^ Np. Consequently, the argument of the Bessel function in Eq. (1T8|) is 
much less than its order. Then HO] 



h{x) 



{x/2f {x/2Ye'' 



r(z/+l) y/2^y^ 



(21) 
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and the maximum probability at low energy is obtained from 



max 



Np/2 



27i{Np + Np 



-1/2 



(22) 



where we use Np ± Np ^ Np. 



Thus, at high as well as at low energies, the maximal value of the net proton 
probability distribution 



1 



27r(iVp + Np 



1/2 



(23) 



is approximately equal and determined by the total mean number of protons 
and antiprotons. 

At intermediate energies, none of the approximations used above hold. How- 
ever, near the maximum, we can approximate P{N) by a Gaussian, 



1 (JV-JV)2 

P{N) ^ -_e-^, (24) 
\'27Ta^ 

with the maximum value ^ 

^"^■^^ ^ 72^' ^^^^ 

In the hadron resonance gas, the variance is directly related to the total num- 
ber of protons and antiprotons 

a^ = Np + Np. (26) 

Consequently, we conclude that the maximum value of the net proton proba- 
bility distribution is, to a good approximation, determined by 

pHRG _ I /^,7^ 

max ~ r _ _ -,1/2 ' J 

27r{Np + Np)' 

irrespective of the collision energy. On the other hand, as P^ax^ depends only 
weakly on the beam energy (see Fig.|3]), the value of Np + Np in the acceptance 
window has to be essentially energy independent. 

r - - -11/2 

In Fig. [5]-right we show the energy dependence of the variance a = Np + Np 
in heavy ion collisions obtained at mid rapidity with \y\ < 0.5, in the trans- 
verse momentum window, 0.4 GeV < pt < 0.8 GeV. The small variation of 
a, seen in this figure, is consistent with the observed weak dependence of the 
maximum of P{N) on A/i^vAf' which thus seems to be a consequence of the 
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Fig. 4. The net-proton probability distribution in central Au-Au collisions at 
200GeV. The HRG results are from Ref. 1351. while the data are due 



to the STAR Collaboration [27]. Also shown are results obtained with the UrQMD 
event generator |41j from a sample of 10^ events. 



particular p^-cut applied to the transverse-momentum distributions of protons 
and antiprotons. 

In the hadron resonance gas, not only the variance, but also all higher moments 
of the net proton distribution can be expressed in terms of the mean number of 
protons and antiprotons. Indeed, using Eqs. ([7]) and (fTOj) . one finds for n = 1, 



X2n 



VT^ 



X2n+1 = [Np - Np) . 



(2J 



Consequently, the skewness S = VT^Xs/ {VT^X2Y^'^ and the kurtosis k 
VT'^Xi/{yT^X2f are directly related to Np and Np by 



S 



Np - Np 
[Np + Npf/^'- 



1 



K 



Np + Np 



(29) 



The energy dependence of skewness and kurtosis is shown in Fig. [5]-right. 
The weak energy dependence of {Np + Np) implies that the kurtosis and all 
even cumulants are essentially independent of beam energy. By contrast, the 
skewness and all odd cumulants are strongly decreasing functions of energy, 
and approach zero at high energies, where Np ~ Np. 

We have seen that in the HRG model all moments of net proton distribution 
can be expressed directly in terms of the mean-number of protons and antipro- 
tons. This can be generalized to moments of the fluctuations of strangeness 
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Fig. 5. Left-hand figure: The mean (M), variance (cr), skewness (S) and kurtosis 
(k) for net proton number in central Au-Au collisions computed using the proba- 
bility distributions shown in Fig. [3l Right-hand figure: The probability distribution 
for net-strangeness in a gas composed of (ir"*", A, H, fi) and their antiparticles at 
T ~ rrij^. The distribution is calculated with and without multi-strange baryons. 



and electric charge. Indeed, using Eqs. (I7j) and ( fTOl) one finds, that in general 



kl 



\g\ 



X2k 



{Nn + Nn)n'\ X2k+ 



2k+l 



(30) 



where 1^1 = 1, 2 or 3 for baryon number, electric charge and strangeness, 
respectively. For the electric charge fluctuations this implies that 



Q_ Ni-N.i + 2{N2-N.2) Q _ iVi + iV-i + 4(iV2 + iV-2) 



^_ Ar^_iV_, + 8(iV,-iV,,) Q 
X3 VT'^ ' 



Ni + iV_i + 16(iV2 + 
VT^ 



, (32) 



where Ni {N-i) and N2 {N_2) are the mean numbers of particles with charge 
+1(— 1) and +2(— 2) respectively. Thus, in the HRG model, the experimental 
input on the mean number of charged and multi-charged particles is suffi- 
cient to compute the net charge probability distribution and the correspond- 
ing moments. We stress, however, that this is only the case in the Boltzmann 
approximation, which in the case of electric charge is a poor approximation. 
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4 The influence of multi-charged particles and quantum statistics 

on P{N) 

The influence of multi-charged particles on fluctuations was discussed in Ref. 

In the HRG model, the product of kurtosis and variance of the net baryon 
number, ko"^, which receives contributions from singly charged particles only, 
is always equal to unity. This result does not depend on the hadron mass 
spectrum nor on the thermal parameters. However, for fluctuations of the 
electric charge, the contribution of particles with charge two, and A , 
leads to a non-trivial T and fi dependence of this quantity [TT] . Furthermore, a 
straightforward calculation shows that in an ideal pion gas, na"^ is temperature 
dependent and differs from unity, due to quantum statistics effects. 

Clearly these effects should be reflected in the corresponding probability dis- 
tribution. To illustrate the influence of multi-charged particles and of quantum 
statistics we consider two simple models: i) a. gas composed of A, S and fl and 
their antiparticles and ii) an ideal pion gas. Within the flrst model we com- 
pute P{S) and within the second P{Q). In Figs. O-right and|B]-right we show 
the probability distributions of the net strangeness and net electric charge 
calculated in the models deflned above. 

To assess the influence of multi-charged particles, we compare the results with 
P{N) computed using Eq. f fTSj) . which accounts only for singly charged par- 
ticles. Since the multi-charged particles are heavy, their contribution to the 
mean-charge is negligible. Consequently, the position of the maximum of P{S) 
is not changed when these states are included. However, the contribution of 
double- and triple-charged particles modifles the tail of the probability dis- 
tribution considerably, and therefore also the higher order moments. Thus, 
a comparison of fluctuations in the hadron resonance gas model with data, 
must account for multi-charged particles, since they are included in the data 
through their decay products. 

The probability distributions in Eqs. f lT6] - [T8|l are valid only in Boltzmann 
approximation. This approximation is justifled for the net strange and net 
baryon number probability distributions. However, owing to the small mass of 
the pion, quantum statistics must be included in studies of charge fluctuations. 
The effect of the corrections to the Boltzmann approximation due to quantum 
statistics is similar to that induced by multi-charged states, and could affect 
the net charge distribution P{Q). 

Indeed, the pressure of charged pions at temperature T and electric charge 
chemical potential /i can be expanded in powers of the fugacity, e^^/"^ ^ 

P7r++P7r = di{f3miYK2{f3mi) cosh(g//3/i), (33) 

^ 1=1 
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Fig. 6. Left-hand figure: The ratios and Rq from Eq. (137p in the HRG model along 
the freeze-out line. Right-hand figure: The probability distributions for net-electric 
charge in a pion gas at T ~ m-,^ calculated with quantum statistics (Exact) and in 
the Boltzmann approximation. 



where mi = lm.„, di = 1//^ and qi = I. Thus, the pion pressure, Eq. ( l33l) . can 
be viewed as a sum of Boltzmann contributions owing to particles with charge 
qi mass mi and degeneracy factor di. 



At /i = 0, the pion contribution to the n-th order cumulant of the electric 
charge reads 

™2 oo 

X^ = ^,T.^''-'mPm^). (34) 

An approximation to Xn obtained by truncating the series flM|) at some 
order 1 = 1*. The value of /*, which yields a good approximation to the exact 
result, depends strongly on the order of the cumulant. 

At temperatures T ~ the dominant contribution to the cumulants Xn 
with n < 5 is due to the / = 1 and / = 2 terms. In other words, in the HRG 
model, the next to leading order Boltzmann approximation is sufficient to 
quantify cumulants up to fifth order. Consequently, Eqs. and (152]) remain 
valid if the effective charge-2 pion contribution in the expansion flMl) is taken 
into account. However, this is not the case for the sixth- and higher-order 
cumulants. The dominant contribution to is due to the / = 3 term in 
Eq. ( I34l) . Consequently, Eq. ( I30l) implies that the sixth-order cumulant of the 
electric charge is well approximated by 



VT^ 



N^ + Ni + 2\N2 + N2. + N-^ + N2^) + t{N^^ + N^^) , (35) 



where A''i (A^-i) and (^-2) are the numbers of particles with charge +1(— 1) 
and +2(— 2) respectively, while Nk-K and iVfc^ are the /c-th order corrections to 
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the Boltzmann approximation for the charged pions numbers 



The number of antipions, N^^, is obtained form Eq. (136|) by the replacement 
/i — — /i, where is the electric charge chemical potential. 

The influence of quantum statistics on fluctuations of the electric charge can 
be made more transparent by considering the ratios: 

r)Q A.5 A.1 pQ A6 A.2 ('3 7\ 

-"-5 — ~~Q Q ' -"-6 — ~Q Q- J 

If on/?/ doubly charged contributions are included in Eqs. (I301 - I321) and (l35l) . the 
ratios -R^ = -R^ = 5 are independent of temperature and chemical potential. 
Any deviation from this value signal contributions emanating from the third- 
and higher-order corrections owing to quantum statistics. 

In Fig. EJ-left, the effect of quantum statistics is illustrated by the ratios -R5 
and Rq obtained in the HRG model along the chemical freeze-out line. As 
shown, the ratio is indeed within < 1 % consistent with the result from 
Eqs. f p0] - [32|) . while Rq deviates by up to 20 % due to the missing effective 
triply charged pionic contribution to Xg ■ 

The effect of quantum statistics is also reflected in the net electric charge 
probability distribution. In Fig. [6]-right we show that the inclusion of quantum 
statistics results in a broadening of the net charge probability distribution of 
an ideal pion gas, which is indeed very similar to the effect of multi-charged 
particles. 

In the HRG the probability distribution of the electric charge, due to the 
relevant contributions from the second and the third order pionic quantum 
statistics corrections, can be explicitly calculated by applying Eq. f|T6|) instead 
of Eq. (P 



5 Conclusions 



We have analyzed the properties of the net charge probability distributions 
in heavy-ion collisions within the hadron resonance gas model. Within this 
model, these distributions and the corresponding cumulants are linked directly 
to other observables, and can be determined from the mean number of charged 
and multi-charged particles. 
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Differences between the probability distributions for strangeness, electric charge 
and net-baryon number were discussed and the role of multi-charged particles 
and quantum statistics was explored. A comparison of the UrQMD results for 
the net baryon probabihty distribution in central Au-Au collisions at x/s^vjv ~ 
200 GeV with data and with the hadron resonance gas model was also pre- 
sented. We find that the distribution obtained within the UrQMD model is 
much narrower than both the experimental distribution and that of the hadron 
resonance gas. 

The comparison of model results with data on net proton distributions and 
the corresponding moments is particularly transparent. This is because, within 
the hadron resonance gas model, these quantities are determined entirely by 
the measured yields of protons and antiprotons. This offers an unambiguous 
route for confronting this particular model with experimental data. 

We have presented model predictions for net proton probability distributions 
P{N) in heavy ion collisions at RHIC and at LHC energies and computed the 
centrality dependence of P{N) in Au-Au collisions at y^^vAf ~ 62.4 and 130 
GeV. The results for P{N) obtained in the hadron resonance gas along with 
the chemical freeze-out curve can be compared with data. Such a comparison 
would be an important test of the relation between the chemical freeze-out 
and the QCD cross-over line in heavy ion experiments. 
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